Gray whale (Eschrichtius robustus) post-mortem findings from December 2018 through 2021 during the Unusual Mortality Event in the Eastern North Pacific.

Beginning in December 2018, increased numbers of gray whale (Eschrichtius robustus) strandings were reported along the west coast of Mexico, the United States, and Canada, prompting declaration of a gray whale Unusual Mortality Event (UME) by the United States National Marine Fisheries Service. Although strandings declined in 2020 and 2021 from a peak in 2019, the UME is still ongoing as of fall 2023. Between 17 December 2018 and 31 December 2021, 503 animals stranded along the west coast of North America, with 226 strandings in Mexico, 71 in California, 12 in Oregon, 56 in Washington, 21 in British Columbia, and 117 in Alaska. These included 187 males, 167 females, and 149 whales of undetermined sex; and 193 adults, 194 subadults, 40 calves, 1 fetus, and 75 whales of undetermined age class. We report on 61 of the 503 carcasses (12%) that had external and internal gross necropsy and/or histopathology data: of these 61 whales, findings that contributed to death were identified in 33 (54%) whales. Sixteen of the 61 (26%) were severely emaciated. Gross lesions of blunt force trauma consistent with vessel strike were identified in 11 of the 61 animals (18%), only two of which were emaciated. Two whales (3%) were entangled at time of death, and one died from entrapment. Signs of killer whale (Orcinus orca) interaction were documented in 19 of the 61 animals; five were deemed from recent interactions and three (5%) likely contributed to mortality. A specific cause of death could not be identified in 28 of 61 whales (46%). Additionally, logistical challenges and the advanced state of decomposition of most examined carcasses precluded detection of potential infectious or toxic causes of morbidity or mortality. Up to 2016, the eastern North Pacific population of gray whale population had generally been increasing since the cessation of historic whaling and a prior UME in 1999-2000. However, recent abundance and calf production estimates have declined, a trend that overlaps the current UME. The relative contributions of carrying capacity, environmental change, prey shifts, and infectious, toxic, and other processes to the increased gray whale mortalities have not yet been resolved. Nevertheless, the marked temporal increase in strandings, including findings of malnutrition in some of the whales, along with low calf production, likely represent consequences of complex and dynamic ecological interactions in the ocean impacting the population.

Bellwethers of change: population modelling of North Pacific humpback whales from 2002 through 2021 reveals shift from recovery to climate response.

For the 40 years after the end of commercial whaling in 1976, humpback whale populations in the North Pacific Ocean exhibited a prolonged period of recovery. Using mark-recapture methods on the largest individual photo-identification dataset ever assembled for a cetacean, we estimated annual ocean-basin-wide abundance for the species from 2002 through 2021. Trends in annual estimates describe strong post-whaling era population recovery from 16 875 (± 5955) in 2002 to a peak abundance estimate of 33 488 (± 4455) in 2012. An apparent 20% decline from 2012 to 2021, 33 488 (± 4455) to 26 662 (± 4192), suggests the population abruptly reached carrying capacity due to loss of prey resources. This was particularly evident for humpback whales wintering in Hawai'i, where, by 2021, estimated abundance had declined by 34% from a peak in 2013, down to abundance levels previously seen in 2006, and contrasted to an absence of decline in Mainland Mexico breeding humpbacks. The strongest marine heatwave recorded globally to date during the 2014-2016 period appeared to have altered the course of species recovery, with enduring effects. Extending this time series will allow humpback whales to serve as an indicator species for the ecosystem in the face of a changing climate.

Behavioural responses of fin whales to military mid-frequency active sonar.

The effect of active sonars on marine mammal behaviour is a topic of considerable interest and scientific investigation. Some whales, including the largest species (blue whales, Balaenoptera musculus), can be impacted by mid-frequency (1-10 kHz) military sonars. Here we apply complementary experimental methods to provide the first experimentally controlled measurements of behavioural responses to military sonar and similar stimuli for a related endangered species, fin whales (Balaenoptera physalus). Analytical methods include: (i) principal component analysis paired with generalized additive mixed models; (ii) hidden Markov models; and (iii) structured expert elicitation using response severity metrics. These approaches provide complementary perspectives on the nature of potential changes within and across individuals. Behavioural changes were detected in five of 15 whales during controlled exposure experiments using mid-frequency active sonar or pseudorandom noise of similar frequency, duration and source and received level. No changes were detected during six control (no noise) sequences. Overall responses were more limited in occurrence, severity and duration than in blue whales and were less dependent upon contextual aspects of exposure and more contingent upon exposure received level. Quantifying the factors influencing marine mammal responses to sonar is critical in assessing and mitigating future impacts.

Boom-bust cycles in gray whales associated with dynamic and changing Arctic conditions.

Climate change is affecting a wide range of global systems, with polar ecosystems experiencing the most rapid change. Although climate impacts affect lower-trophic-level and short-lived species most directly, it is less clear how long-lived and mobile species will respond to rapid polar warming because they may have the short-term ability to accommodate ecological disruptions while adapting to new conditions. We found that the population dynamics of an iconic and highly mobile polar-associated species are tightly coupled to Arctic prey availability and access to feeding areas. When low prey biomass coincided with high ice cover, gray whales experienced major mortality events, each reducing the population by 15 to 25%. This suggests that even mobile, long-lived species are sensitive to dynamic and changing conditions as the Arctic warms.

Migratory destinations and spatial structuring of humpback whales (Megaptera novaeangliae) wintering off Nicaragua.

Understanding the migratory patterns of large whales is of conservation importance, especially in identifying threats to specific populations. Migration ecology, including migratory destinations, movements and site fidelity for humpback whales (Megaptera novaeangliae) remain poorly studied in parts of the range of the Central America population, considered endangered under the United States Endangered Species Act. This study aimed to investigate the migratory destinations of humpback whales sighted at two study sites in Nicaragua, which are part of the Central America population. A ten-year photographic database of humpback whales observed off Nicaragua was combined with citizen science contributions and sightings from dedicated research programs. The resulting image collection was compared with available historical photo identifications and databases using an automated image recognition algorithm. This approach yielded 36 years of photographic identification totaling 431 recaptures in Nicaragua (2006-2008 and 2016-2021) and 2539 recaptures (1986-2020) in both feeding and breeding grounds of 176 unique individuals sighted in Nicaragua. Our results showed that photo-identified whales were recaptured between October and April in breeding grounds and year-round in feeding grounds between British Columbia and California, with peak recaptures between June and October. Our study provided first-time evidence on fine-scale site affinity of individual humpback whales within Nicaraguan waters and to other breeding and feeding grounds.

A collaborative and near-comprehensive North Pacific humpback whale photo-ID dataset.

We present an ocean-basin-scale dataset that includes tail fluke photographic identification (photo-ID) and encounter data for most living individual humpback whales (Megaptera novaeangliae) in the North Pacific Ocean. The dataset was built through a broad collaboration combining 39 separate curated photo-ID catalogs, supplemented with community science data. Data from throughout the North Pacific were aggregated into 13 regions, including six breeding regions, six feeding regions, and one migratory corridor. All images were compared with minimal pre-processing using a recently developed image recognition algorithm based on machine learning through artificial intelligence; this system is capable of rapidly detecting matches between individuals with an estimated 97-99% accuracy. For the 2001-2021 study period, a total of 27,956 unique individuals were documented in 157,350 encounters. Each individual was encountered, on average, in 5.6 sampling periods (i.e., breeding and feeding seasons), with an annual average of 87% of whales encountered in more than one season. The combined dataset and image recognition tool represents a living and accessible resource for collaborative, basin-wide studies of a keystone marine mammal in a time of rapid ecological change.

Minke whale feeding rate limitations suggest constraints on the minimum body size for engulfment filtration feeding.

Bulk filter feeding has enabled gigantism throughout evolutionary history. The largest animals, extant rorqual whales, utilize intermittent engulfment filtration feeding (lunge feeding), which increases in efficiency with body size, enabling their gigantism. The smallest extant rorquals (7-10 m minke whales), however, still exhibit short-term foraging efficiencies several times greater than smaller non-filter-feeding cetaceans, raising the question of why smaller animals do not utilize this foraging modality. We collected 437 h of bio-logging data from 23 Antarctic minke whales (Balaenoptera bonaerensis) to test the relationship of feeding rates (λf) to body size. Here, we show that while ultra-high nighttime λf (mean ± s.d.: 165 ± 40 lunges h-1; max: 236 lunges h-1; mean depth: 28 ± 46 m) were indistinguishable from predictions from observations of larger species, daytime λf (mean depth: 72 ± 72 m) were only 25-40% of predicted rates. Both λf were near the maxima allowed by calculated biomechanical, physiological and environmental constraints, but these temporal constraints meant that maximum λf was below the expected λf for animals smaller than ~5 m-the length of weaned minke whales. Our findings suggest that minimum size for specific filter-feeding body plans may relate broadly to temporal restrictions on filtration rate and have implications for the evolution of filter feeding.

Influence of Life-History Parameters on Persistent Organic Pollutant Concentrations in Blubber of Eastern North Pacific Gray Whales (Eschrichtius robustus).

Exposure to persistent organic pollutants (POPs) can significantly impact marine mammal health, reproduction, and fitness. This study addresses a significant 20-year gap in gray whale contaminant monitoring through analysis of POPs in 120 blubber biopsies. The scope of this substantial sample set is noteworthy in its range and diversity with collection between 2003 and 2017 along North America's west coast and across diverse sex, age, and reproductive parameters, including paired mothers and calves. Mean blubber concentrations of polychlorinated biphenyls (∑PCBs), dichlorodiphenyltrichloroethanes (∑DDTs), and chlordanes (∑CHLs) generally decreased since previous reports (1968-1999). This is the first report of polybrominated diphenyl ethers (PBDEs) and select hexachlorocyclohexanes (HCHs) in this species. Statistical modeling of the 19 most frequently detected compounds in this dataset revealed sex-, age-, and reproductive status-related patterns, predominantly attributed to maternal offloading. Mean POP concentrations differed significantly by sex in adults (17 compounds, up to 3-fold higher in males) but not in immatures (all 19 compounds). Mean POP concentrations were significantly greater in adults versus immatures in both males (17 compounds, up to 12-fold) and females (13 compounds, up to 3-fold). POP concentrations were detected with compound-specific patterns in nursing calves, confirming maternal offloading for the first time in this species.

Oceanic giants dance to atmospheric rhythms: Ephemeral wind-driven resource tracking by blue whales.

Trophic transfer of energy through marine food webs is strongly influenced by prey aggregation and its exploitation by predators. Rapid aggregation of some marine fish and crustacean forage species during wind-driven coastal upwelling has recently been discovered, motivating the hypothesis that predators of these forage species track the upwelling circulation in which prey aggregation occurs. We examine this hypothesis in the central California Current Ecosystem using integrative observations of upwelling dynamics, forage species' aggregation, and blue whale movement. Directional origins of blue whale calls repeatedly tracked upwelling plume circulation when wind-driven upwelling intensified and aggregation of forage species was heightened. Our findings illustrate a resource tracking strategy by which blue whales may maximize energy gain amid ephemeral foraging opportunities. These findings have implications for the ecology and conservation of diverse predators that are sustained by forage populations whose behaviour is responsive to episodic environmental dynamics.

Blue whales increase feeding rates at fine-scale ocean features.

Marine predators face the challenge of reliably finding prey that is patchily distributed in space and time. Predators make movement decisions at multiple spatial and temporal scales, yet we have a limited understanding of how habitat selection at multiple scales translates into foraging performance. In the ocean, there is mounting evidence that submesoscale (i.e. less than 100 km) processes drive the formation of dense prey patches that should hypothetically provide feeding hot spots and increase predator foraging success. Here, we integrated environmental remote-sensing with high-resolution animal-borne biologging data to evaluate submesoscale surface current features in relation to the habitat selection and foraging performance of blue whales in the California Current System. Our study revealed a consistent functional relationship in which blue whales disproportionately foraged within dynamic aggregative submesoscale features at both the regional and feeding site scales across seasons, regions and years. Moreover, we found that blue whale feeding rates increased in areas with stronger aggregative features, suggesting that these features indicate areas of higher prey density. The use of fine-scale, dynamic features by foraging blue whales underscores the need to take these features into account when designating critical habitat and may help inform strategies to mitigate the impacts of human activities for the species.

Scaling of maneuvering performance in baleen whales: larger whales outperform expectations.

Despite their enormous size, whales make their living as voracious predators. To catch their much smaller, more maneuverable prey, they have developed several unique locomotor strategies that require high energetic input, high mechanical power output and a surprising degree of agility. To better understand how body size affects maneuverability at the largest scale, we used bio-logging data, aerial photogrammetry and a high-throughput approach to quantify the maneuvering performance of seven species of free-swimming baleen whale. We found that as body size increases, absolute maneuvering performance decreases: larger whales use lower accelerations and perform slower pitch-changes, rolls and turns than smaller species. We also found that baleen whales exhibit positive allometry of maneuvering performance: relative to their body size, larger whales use higher accelerations, and perform faster pitch-changes, rolls and certain types of turns than smaller species. However, not all maneuvers were impacted by body size in the same way, and we found that larger whales behaviorally adjust for their decreased agility by using turns that they can perform more effectively. The positive allometry of maneuvering performance suggests that large whales have compensated for their increased body size by evolving more effective control surfaces and by preferentially selecting maneuvers that play to their strengths.

Naturally stressed? Glucocorticoid profiles in blubber of blue and gray whales in response to life history parameters.

The goal of the present study was to carry out a thorough methodological validation and describe baseline profiles for glucocorticoid hormones (cortisol and corticosterone) in blubber from blue (n = 77) and gray (n = 103) whales from the eastern North Pacific Ocean. For each species, we modelled cortisol and corticosterone concentrations in response to life history parameters (age, sex, reproductive status) and season or geographic location. In blue whales, cortisol concentrations did not vary significantly by age class, sex, or reproductive status, whereas corticosterone was significantly lower in immature than in adult females (p < .001). In gray whales, cortisol concentrations were significantly higher in lactating whales (p < .05), while corticosterone was significantly different between females and males (p = .001) and elevated in calves (p = .003). In gray whales, corticosterone concentrations were significantly lower in males sampled later in the year (August to November) compared to both sexes sampled between March and August (p = .05), but no seasonal trend occurred in blue whales. Our results indicate that glucocorticoid actions vary between species and sex in large whales. Analysis of multiple hormones improves our understanding of the physiology of maintaining metabolic homeostasis or coping with chronic stressors.

Baleen whale prey consumption based on high-resolution foraging measurements.

Baleen whales influence their ecosystems through immense prey consumption and nutrient recycling1-3. It is difficult to accurately gauge the magnitude of their current or historic ecosystem role without measuring feeding rates and prey consumed. To date, prey consumption of the largest species has been estimated using metabolic models3-9 based on extrapolations that lack empirical validation. Here, we used tags deployed on seven baleen whale (Mysticeti) species (n = 321 tag deployments) in conjunction with acoustic measurements of prey density to calculate prey consumption at daily to annual scales from the Atlantic, Pacific, and Southern Oceans. Our results suggest that previous studies3-9 have underestimated baleen whale prey consumption by threefold or more in some ecosystems. In the Southern Ocean alone, we calculate that pre-whaling populations of mysticetes annually consumed 430 million tonnes of Antarctic krill (Euphausia superba), twice the current estimated total biomass of E. superba10, and more than twice the global catch of marine fisheries today11. Larger whale populations may have supported higher productivity in large marine regions through enhanced nutrient recycling: our findings suggest mysticetes recycled 1.2 × 104 tonnes iron yr-1 in the Southern Ocean before whaling compared to 1.2 × 103 tonnes iron yr-1 recycled by whales today. The recovery of baleen whales and their nutrient recycling services2,3,7 could augment productivity and restore ecosystem function lost during 20th century whaling12,13.

Epizootiology of a Cryptococcus gattii outbreak in porpoises and dolphins from the Salish Sea.

Cryptococcus gattii is a fungal pathogen that primarily affects the respiratory and nervous systems of humans and other animals. C. gattii emerged in temperate North America in 1999 as a multispecies outbreak of cryptococcosis in British Columbia (Canada) and Washington State and Oregon (USA), affecting humans, domestic animals, and wildlife. Here we describe the C. gattii epizootic in odontocetes. Cases of C. gattii were identified in 42 odontocetes in Washington and British Columbia between 1997 and 2016. Species affected included harbor porpoises Phocoena phocoena (n = 26), Dall's porpoises Phocoenoides dalli (n = 14), and Pacific white-sided dolphins Lagenorhynchus obliquidens (n = 2). The probable index case was identified in an adult male Dall's porpoise in 1997, 2 yr prior to the initial terrestrial outbreak. The spatiotemporal extent of the C. gattii epizootic was defined, and cases in odontocetes were found to be clustered around terrestrial C. gattii hotspots. Case-control analyses with stranded, uninfected odontocetes revealed that risk factors for infection were species (Dall's porpoises), age class (adult animals), and season (winter). This study suggests that mycoses are an emerging source of mortality for odontocetes, and that outbreaks may be associated with anthropogenic environmental disturbance.

Application of endocrine biomarkers to update information on reproductive physiology in gray whale (Eschrichtius robustus).

Most of our knowledge on reproductive biology of gray whales dates back to scientific research conducted during commercial whaling in the late 1950s and 1960s. The goal of the present study was to provide updated insights on reproductive physiology of gray whales, using progesterone and testosterone as biomarkers. We measured hormone concentrations using enzyme immunoassay (EIA) techniques in blubber biopsies collected from 106 individual whales from March to November over a span of 12 years (2004-2016) between California and Alaska. We found testosterone concentrations in males to increase significantly with age (P = 0.03). Adult males showed significantly elevated testosterone concentrations when sampled in the fall compared to the summer (P = 0.01), likely indicating physiological preparation for mating. We measured testosterone concentrations in females of different age classes, but no statistical differences were found. We found significantly higher progesterone concentrations in pregnant females compared to non-pregnant females and adult males (P< 0.001), indicating progesterone is a valid biomarker for pregnancy in gray whales. Both female and male calves had elevated progesterone concentrations, suggesting maternal transfer via lactation. We fit a mixture of two normal distributions to progesterone data from all non-calf females to identify clusters of high and low progesterone and estimated the probability of being pregnant for whales of unknown reproductive status. With this approach we identified likely pregnant and non-pregnant animals. This study represents an important milestone on reproductive profiles in this population, that can be used to estimate more accurate and precise reproductive parameters to be used for better understanding population dynamics of gray whales.

Blubber endocrine profiles provide insights into reproductive biology in blue whales from the eastern North Pacific Ocean.

The goal of the present study was to complement existing data of testosterone and progesterone in blue whale (Balaenoptera musculus) blubber from the eastern North Pacific Ocean to evaluate effects of seasonality and location on these hormones and to better assess reproductive status of individuals. Physiological parameters regarding reproduction are fundamental for describing population dynamics, and hormones can be a valid tool to estimate those for wildlife populations. In this study, blubber tissue was validated for testosterone and progesterone assays. Hormone concentrations were measured in 69 (35 males and 34 females) blubber samples from live (n = 66) and stranded (n = 3) animals collected between 2002 and 2016 from a known winter reproductive ground in the Gulf of California (GoC) and summer feeding areas along the United States West Coast (USWC), specifically off the states of California and Oregon. Results were combined with sighting histories as a tool to determine reproductive status of individual whales. Testosterone concentrations in adult male blue whales were significantly higher (p < 0.05) in blubber biopsies sampled off the USWC between the months of June and November compared to those sampled in the GoC between February and April. Elevated testosterone concentrations likely indicate physiological preparation for reproductive activity while the animals were present off the USWC. Progesterone concentrations were significantly elevated in pregnant females, confirming progesterone as an indicator of pregnancy in blue whales. Probabilities of being pregnant were estimated for adult females with unknown sighting histories based on progesterone concentrations. Testosterone in females was detected and measured only in pregnant whales suggesting its biosynthesis or metabolism is altered during gestation. These results provide updated and new information on the reproductive cycle of blue whales in the eastern North Pacific, posing new milestones to better estimate the timing of the mating season for this endangered population.

Congenital Diseases in Harbor Seals (Phoca vitulina richardsii) from the Salish Sea.

Postmortem data for harbor seals (Phoca vitulina richardsii) in the Salish Sea were analyzed for epidemiologic trends in congenital diseases. Cleft palate, cleft lips, or both (n=8) and cardiac defects (n=5) were the most common congenital abnormalities, followed by cases with multiple defects (n=4). No temporal trends or spatial clusters of cases were seen from 2003 to 2019, during which time monitoring effort was consistent. Cases could not be linked to specific causes such as environmental contamination or maternal malnutrition. Our study suggests that a yearly prevalence of 2.9%±2.2 is the endemic level of congenital disease in this stable harbor seal population. Continued monitoring of birth defects and overall harbor seal population status could help to identify emerging teratogens.

Population structure in a continuously distributed coastal marine species, the harbor porpoise, based on microhaplotypes derived from poor-quality samples.

Harbor porpoise in the North Pacific are found in coastal waters from southern California to Japan, but population structure is poorly known outside of a few local areas. We used multiplexed amplicon sequencing of 292 loci and genotyped clusters of single nucleotide polymoirphisms as microhaplotypes (N = 271 samples) in addition to mitochondrial (mtDNA) sequence data (N = 413 samples) to examine the genetic structure from samples collected along the Pacific coast and inland waterways from California to southern British Columbia. We confirmed an overall pattern of strong isolation-by-distance, suggesting that individual dispersal is restricted. We also found evidence of regions where genetic differences are larger than expected based on geographical distance alone, implying current or historical barriers to gene flow. In particular, the southernmost population in California is genetically distinct (FST  = 0.02 [microhaplotypes]; 0.31 [mtDNA]), with both reduced genetic variability and high frequency of an otherwise rare mtDNA haplotype. At the northern end of our study range, we found significant genetic differentiation of samples from the Strait of Georgia, previously identified as a potential biogeographical boundary or secondary contact zone between harbor porpoise populations. Association of microhaplotypes with remotely sensed environmental variables indicated potential local adaptation, especially at the southern end of the species' range. These results inform conservation and management for this nearshore species, illustrate the value of genomic methods for detecting patterns of genetic structure within a continuously distributed marine species, and highlight the power of microhaplotype genotyping for detecting genetic structure in harbor porpoises despite reliance on poor-quality samples.

Context-dependent variability in the predicted daily energetic costs of disturbance for blue whales.

Assessing the long-term consequences of sub-lethal anthropogenic disturbance on wildlife populations requires integrating data on fine-scale individual behavior and physiology into spatially and temporally broader, population-level inference. A typical behavioral response to disturbance is the cessation of foraging, which can be translated into a common metric of energetic cost. However, this necessitates detailed empirical information on baseline movements, activity budgets, feeding rates and energy intake, as well as the probability of an individual responding to the disturbance-inducing stressor within different exposure contexts. Here, we integrated data from blue whales (Balaenoptera musculus) experimentally exposed to military active sonar signals with fine-scale measurements of baseline behavior over multiple days or weeks obtained from accelerometry loggers, telemetry tracking and prey sampling. Specifically, we developed daily simulations of movement, feeding behavior and exposure to localized sonar events of increasing duration and intensity and predicted the effects of this disturbance source on the daily energy intake of an individual. Activity budgets and movements were highly variable in space and time and among individuals, resulting in large variability in predicted energetic intake and costs. In half of our simulations, an individual's energy intake was unaffected by the simulated source. However, some individuals lost their entire daily energy intake under brief or weak exposure scenarios. Given this large variation, population-level models will have to assess the consequences of the entire distribution of energetic costs, rather than only consider single summary statistics. The shape of the exposure-response functions also strongly influenced predictions, reinforcing the need for contextually explicit experiments and improved mechanistic understanding of the processes driving behavioral and physiological responses to disturbance. This study presents a robust approach for integrating different types of empirical information to assess the effects of disturbance at spatio-temporal and ecological scales that are relevant to management and conservation.

Isotope-based inferences of the seasonal foraging and migratory strategies of blue whales in the eastern Pacific Ocean.

Migratory marine megafauna generally move vast distances between productive foraging grounds and environmentally stable breeding grounds, but characterizing how they use these habitats to maintain homeostasis and reproduce is difficult. We used isotope analysis of blue whale skin strata (n = 621) and potential prey (n = 300) to examine their migratory and foraging strategies in the eastern Pacific Ocean. Our results suggest that most whales in the northeast Pacific use a mixed income and capital breeding strategy, and use the California Current Ecosystem as their primary summer-fall foraging ground. A subset of individuals exhibited migratory plasticity and spend most of the year in the Gulf of California or Costa Rica Dome, two regions believed to be their primary winter-spring breeding grounds. Isotope data also revealed that whales in the southern Eastern Tropical Pacific generally do not forage in the northeast Pacific, which suggests a north-south population structure with a boundary near the equator.